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Abstract 

We demonstrate that the free motion of any two-dimensional rigid body col¬ 
liding elastically with two parallel, flat walls is equivalent to a billiard system. 
Using this equivalence, we analyze the integrable and chaotic properties of 
this new class of billiards. This provides a demonstration that coin tossing, 
the prototypical example of an independent random process, is a completely 
chaotic (Bernoulli) problem. The related question of which billiard geometries 

can be represented as rigid body systems is examined. 
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I. INTRODUCTION 


Classical dynamics can describe chaotic motion. It is generally accepted that trnly 
chaotic systems are isomorphic to Bernonlli systems. This means that any chaotic solntion 
is eqnivalent in a well-dehned sense to a stochastic one generated by an independent random 
process. The prototype of snch a process is coin tossing described by repeated random trials 
with only two possible ontcomes for each trial with probabilities which do not change dnring 
the trials. These trials are called Bernonlli trials, hence the name. However, the tossing of a 
coin can also be viewed as a simple dynamical system since it is merely the motion of a rigid 
body snbjected to simple forces and bonndary conditions. It is the legitimate to ask whether 
coin tossing is de facto a Bernonlli system when considered as a dynamical process. The 
experience of gamblers throngh the ages wonld snggest that it is. The aim of the present 
paper is to provide a more mathematical demonstration. 

To be specihc, we shall consider a rigid body in two dimensions which moves in the 
absence of forces and makes elastic collisions with two inhnitely massive parallel walls. 
The special case in which this rigid body is simply a stick of zero width will serve to 
model coin tossing. We shall demonstrate that this system is eqnivalent to that of a point 
particle colliding with a snitably cnrved wall. The key featnre to this eqnivalence is the 
fact that elastic collisions by the rigid body necessarily imply specular reflections for the 
related point particle, i.e., the behavionr of a traditional billiard. For the special case of 
tossing a coin between plane walls, the eqnivalent billiard conhgnration is that of a particle 
bonncing between two convex sinusoidal boundary walls. It has been shown in |l| that 
billiards with strictly convex walls are of the Bernoulli type indicating that coin tossing, 
when considered as a dynamical system, is truly a completely random process. However, 
for billiards arising from coin tossing, the correlation between the initial orientation of the 
coin and its orientation on subsequent bounces dies out with an exponential envelope. This 
then provides a characteristic decorrelation time after which the process is random. A fair 
coin toss must necessarily allow the coin to bounce repeatedly on the floor and the ceiling. 
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Evidently, a simple toss and catch process, the colloquial example of a random process, does 
not suffice. 

We shall see that every such rigid body system is equivalent to a unique billiard problem. 
The generality of this equivalence is of some interest. For some time, billiard problems 
have been known to provide convenient and illustrative examples of the general properties 
of Hamiltonian dynamical systems [@-@]. The quantum properties of billiards have also 
been studied (see, for instance [||). Physical realizations of billiard systems have been 
found in microwave cavities as well as in electronic nanometric semiconductor devices and 
used to investigate how the transition from regular to chaotic motion influences the wave 
properties of these systems p-§|. Yet, many of the billiard geometries studied are somewhat 
arbitrary and do not represent real physical systems. The present results offer a broad class 
of physically realizable dynamical systems exhibiting all types of motion; integrable, near 
integrable (KAM), and chaotic motion of increasing randomness (iP-flows, C-flows, and 
Bernoulli flows). 

The organization of the paper is as follows. In section |^, we consider the specihc problem 
of the stick (or coin). The general transformations relating the rigid body problem to the 
billiard system are derived, and the condition for specular reflection is proved. Section H 
is deals with some of the properties associated with elliptical rigid bodies. In particular, we 
consider how the onset of chaos depends on the eccentricity of the ellipse and the separation 
between the walls. In section we consider the inverse problem of deriving body shapes 
from given wall shapes. Finally, a number of conclusions are drawn in section ^ along 
with some suggestions for certain interesting shapes and possible directions for a quantum 
mechanical analysis of this problem. 

II. TRANSFORMATION EQUATIONS 
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A. Sticks 


Let us first consider the simple case of a stick of length 2L, mass M, and radius of gyration 
K, bouncing freely between two flat walls separated by a distance H. (For simplicity, we 
assume that the center of mass is at the mid-point of the stick.) The position of the stick is 
characterized by y, the height of the center of mass above the lower wall, and 6, the angle 
of rotation of the stick from the vertical. If we introduce the dimensionless height y = y/ k, 
the (scaled) energy of the stick is given by 

B = Y(9"+y). (1) 

This looks formally like the energy of a free point particle moving in an Euclidean plane 
parameterized by the dimensionless coordinates Xi = 9 and X 2 = rj- At the collision of the 
stick with the wall, it is easy to see that 

{MK)Sr] = fnSt, 

{Mk^) 69 = Lsm{6)fnSt, (2) 


where /„ is the normal force at the wall. This gives a linear relation between the change in 
velocities at the collision point. 


se 


L 

R 


sin{d)Sri. 


( 3 ) 


Now, the boundary for the motion of our point particle is given by the minimum value of 
f] for a given angular orientation. From conservation of energy we hnd (using (1) and (3)) 
that 


. ^ -2{r] + {L/K)esm{e)) 

^ i + {L/Kysm{ey ’ 

resulting in r^min = {L/k,) \ cos( 6')|. Therefore, the bottom and top boundary walls {b{9) and 
t{6), respectively) corresponding to the billiard problem for the stick are 

b(9) = —I cos(6')|, 

K 

t{9) = H -|cos(0)|. (5) 

K 
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Finally, consider the infinitesimal change in the velocity vector of our point particle after 
a collision, 

6v = {f-6v)f+{N-6v)N. (6) 

Here, T and N are the orthonormal tangent and normal vectors at the billiard boundary. 
The tangent vector along the lower wall is 

f =(l,--sm(«)). (7) 

K 

Using (3), we have 

f ■ (5F = f ■ (50, 5?)) = 0, (8) 

and consequently, the tangential component of 5v vanishes. Since the magnitude of v is 
preserved by energy conservation, the above result indicates that upon collision, the normal 
component of 5v simply changes sign. Thus, the point particle in our equivalent undergoes 
specular reflections at the boundary walls. Hence, the motion of a rigid stick is equivalent 
to a billiard problem with properly chosen boundaries. (Since the horizontal 0-axis is clearly 
periodic, one may consider the configuration space to be the surface of a cylinder with 
‘undulating’ top and bottom edges.) 

B. Arbitrary Shapes and Mass Distributions 

It is a simple matter to generalize the above to the case of rigid bodies of arbitrary 
shape. We restrict our discussion to objects possessing at least one symmetry axis in the 
plane. (The procedure for asymmetric objects is identical except that the upper boundary 
is shifted by tt with respect to the lower.) Consider a closed plane curve C. If C is not 
everywhere convex, we replace C with its convex envelope since concave segments cannot 
come into contact with the boundary walls. As before, the inertial properties of the object 
will be expressed in terms of the total mass M and the radius of gyration k. The height of 
the center of mass of the object above the wall will be denoted by rj (scaled by the gyration 
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factor K,). The curve C can be specified completely by its radius, R{a), as a function of the 
angle a from some body-hxed axis d. 


R{a) = (-R(q;) cos(q;), -R(q;) sin(Q;)) 


It is elementary to construct the tangent to C, dehned as 


(9) 


T(q;) = dR/da , 


( 10 ) 


and the associated normal vector N{a). When the point of contact between the body and 
the wall is at ac, the perpendicular distance between the center of mass and the wall is 
simply 

1 


o(ac) = --=- , I 

K {N-NyR « VR^ + R'^ 


( 11 ) 


To obtain the desired wall shape, we must determine that orientation of the rigid body, 6, 
which makes ac the point of contact. (For convenience, we measure 6, clockwise, from the 
normal to the bottom wall to the body-fixed axis a.) Obviously, the condition of tangency 
is that N{ac) is parallel to the wall normal: 

R'(ac) 


R{ac 


= tan(Q;c + 6^), 


( 12 ) 


which then produces the lower boundary h{d) = r]{ac{0)) with the upper wall being given 
by t{9) = H — b{6) (or t{6) = H — b{6) + n for asymmetric objects). 

The equations describing the actual collision have a form similar to equations (2), 


= fnSt, 

{MK^dO = R{ac) sin(Q;c -|- 6)fnSt, 


(13) 


with the energy again being given by (1). Since equation (13) follows from (2) with the 
substitutions 


L R{ac) , 

9^9 + ac , (14) 
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it is clear that the proof for specular reflection is identical to that given for the stick. Thus, 
for every shape of the rigid body, there is an equivalent billiard problem. 

As an example of the above procedure, consider an elliptical disk with a semi-major 
axis of one and a semi-minor axis of l/(Vl -f e). Imagine that the center of mass is at the 
geometric center of the ellipse. Relative to this point, the ellipse is parameterized as 


R{a) = 


(1 -I- esin^(Q;))^A 


Equation (12) leads to 


cos(q;c + 9) = 


(1 -|- esin^(Q;c)) 

(1 -|- (e^ -|- 2e) sin^(Q;c))^/^ 


and 


(15) 


(16) 


sm ctr = 


sin2(0) 


(17) 


(1 -f- (e^ -|- 2e) 008^(9)) 

Substituting (15), (16), and (17) into (11) and (12) results in the following billiard boundary 

' 1 \ /(I -|- e cos'^{9)) 


m = 


K 


1 + 6 


(18) 


In the limit of e —0, h{9) 1 /k which is the obvious result for a circle of radius one. In the 

other limit e ^ oo, the ellipse degenerates into a stick and we recover the previous result of 
b{9) = {1 /k) \ cos( 


III. ELLIPTIC DISCS 

KAM motion: We now focus on the special case in which the rigid body is an ellipse so 
that the equivalent billiard problem has the wall described by equation (18) of the preceding 
section. A Poincare section map for the elliptic system with e = 0.003 and 77 = 3 is shown 
in Fig.2. This is a phase portrait of the contact point 9c versus the angle of incidence <f)c 
with respect to the local normal at 9c- The regular and chaotic regions of this “mixed” 
system are clearly visible. For such a small value of the eccentricity e, the boundary can be 
expanded as 
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6(»)~(l-(£/2)sm^(»)). 


(19) 


This indicates that even a small variation from the integrable billiard generated by a circle 
can produce a dense population of stochastic trajectories. The corresponding phase space 
can be described using the KAM theorem [Q. For instance, Fig.3 is a phase portrait for 
e = 0.55 and H = 1.9. This section map contains both hyperbolic and elliptical regions 
of phase space. The inner circles correspond to KAM tori of irrational winding number 
frequencies. The island chains encircling these tori surround the elliptic fixed points of the 
motion. As is well known, this structure is self-similar and persists at smaller length scales 
ad infinitum. 

The transition from a mixed geometry to a chaotic one occurs suddenly in these rigid 
body problems. As long as the motion is purely chaotic, the actual value of H plays no 
explicit role. When the motion contains both hyperbolic and elliptic regions, the relative 
size of these different regions is dictated by the height H with the elliptic regions disappearing 
as H and/or e become larger. It is somewhat delicate to make a numerical determination of 
the curve which separates the KAM and chaotic regions. This is particularly true when H 
is large and e is small. The approximate boundary curve for the elliptic problem is 
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17 

IP 


( 20 ) 


This expression reproduces the results of simulations to about 1% which is probably the limit 
of their accuracy. In the limit as e —> 0 and the object approaches a circle, the wall function 
approaches a constant plus a small sinusoidal oscillation. Given any e, equation (20) shows 
us that it is always possible to pick H sufficiently large that K-type motion results. 

Application to coins: The case of a very eccentric ellipse can be viewed as a model of a 
coin tossing process. Regarding this ellipse as a coin observed edgewise, we place an oriented 
arrow along its length. With this fixed orientation, we assign the condition “heads” when 
the arrow points to the right and “tails” when the arrow points to the left. The motion of the 
coin consists of the one-dimensional vertical motion of the center of mass and a superimposed 
rotation. This leads to a strictly convex billiard problem b{6) = |cos(6*)| which has been 
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shown to be isomorphic to a Bernoulli shift. Surprisingly, this includes the case where 
the distance between the walls is less than the length of the stick. (In our coin analogy, this 
means that the room is not high enough to permit the coin to flip from heads to tails.) In 
this case, the full angular region is not accessible. Nevertheless, the motion is chaotic. This 
is consistent with equation (20) which indicates that, for any H, it is always possible to pick 
e sufficiently large that K-type motion results. 

In spite of its common use as a random process, the traditional tossing of a coin (without 
bounces) is nothing of the sort. As we have shown, the present coin toss (including collisions 
with a wall) does result in K-type motion. In this regard, it is interesting to consider the 
probability, P{n), that, if the coin shows a “head” at a given collision with the bottom wall, it 
will also show a “head” n collisions later. We expect P{n) —1/2 for an independent random 
process. We have studied P{n) numerically as a function of P[ for a very eccentric ellipse 
(e = 99). It is readily seen that P{n) approaches the asymptotic value of 1/2 exponentially. 

F(n) - 1/2 ~ exp {-n/T{H)) (21) 

When the distance between the walls is equal to the diameter of the coin, the correlation 
time r is evidently inhnite. As the wall separation is increased slightly, the correlation time 
is approximately inversely proportional to the difference between H and the diameter of 
the coin. For example, when H is 10% larger than the diameter of the coin, 20 bounces 
are sufficient to reduce the non-asymptotic part of P{n) to 0.01. When p[ is 20% larger 
than the diameter, 12 bounces suffice. It requires exceptionally good statistics to follow the 
non-asymptotic part of P{n) when p[ becomes appreciable. The “house rules” at casinos 
inevitably require that dice be bounced from a wall. The present results suggest that this 
bounce plays a signihcant role in randomizing the outcome and is not merely a time-honored 
tradition. 


9 


IV. FROM WALL SHAPES TO BODY SHAPES 


We have demonstrated that every problem involving elastic collisions between a rigid 
body and parallel walls is equivalent to a billiard with a suitable (and unique) choice of the 
wall function b{6). Here, we wish to address the inverse problem. For convenience, we shall 
restrict our attention to periodic walls which are everywhere differentiable. (This is not to 
deny that rigid bodies with sharp corners are also of interest.) There are several questions 
of interest: How does one determine R{a) given b{d)l Does every periodic b{9) correspond 
to a realizable rigid body and, if so, is its shape unique? The general technique for going 
from b{9) to R{a) is easy to state but can be challenging to implement. Formally, we use 
equation (12) to eliminate 9 from equation (11). This leads to a non-linear but first-order 
differential equation for R{a). Two observations are useful. First, the extrema in b (as a 
function of 9) coincide with the extrema of R (as a function of a) with a and 9 related by 
equation (12). Second, at such extrema, a = 9 and R{a) = b{9). While it may be difficult 
to hnd analytic solutions to the resulting equation, numerical solutions are accessible. It is 
not always easy to hnd them. For example, imagine that the wall shape is given by equation 
(17). The procedure stated leads to a quadratic equation for R'{a) which is readily solved 
to yield 

R' —e sin a cos a ± 1 — e-|-(l-|-e)(l-|-e sin^ a)R^ 

~D ~ i i 

R 1 -|- e sm q; 

which is consistent with equation (15). Note, however, this consistency is realized by two 
degenerate solutions which renders the numerical problem delicate. 

As we shall see, every periodic b{9) can be regarded as an equivalent rigid body problem. 
According to a well-known theorem, b{9) must have at least four extrema to represent a 
physical wall shape. Every non-constant, periodic function must have at least one maximum 
and one minimum. Thus, the four-extremum requirement can always be met by mapping 
9 9/2. Before proceeding, we note a fundamental ambiguity in the inverse problem. 

If one is given only the equivalent billiard problem, the distinction between b{9) and FT, 
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the separation between the parallel walls, cannot be made uniquely. Thus, b{d) can only be 
determined up to some constant shift Bq. Due to the non-linear nature of the inverse problem, 
the corresponding R{a) generally depends on the value of bo assumed. Hence, the equivalent 
billiard does not determine the rigid body uniquely. Moreover, there will generally be choices 
of bo for which no physically realizable R{a) exists. (The specihc example of elliptic discs is 
considered in the Appendix where we show that no physical solution exists when bo is less 
than some critical value.) 

While the freedom of choosing bo and the associated multiplicity of body shapes which 
arise from its exploitation may seem a unfortunate, there is a silver lining. This freedom 
allows us to see that there are always (inhnitely many) rigid body shapes corresponding to 
any choice of b{6). Proof of existence is compensation for loss of uniqueness. Consider the 
elliptic example for large, positive bo so that ^-dependent contributions to the wall function 
are small relative to bo. The angles a and 6 are almost equal, and the resulting body is 
almost a circle. For sufficiently large bo, this will hold for all b{6). Then, 

bo + b (a + tan“^ [ ^ ^ = (23) 

V \^JJ 

with 

R = bo + r{a) . (24) 

It is understood that bo » r. This inequality entitles us to expand equation (28) in powers 
of r/bo and r'/bo which yields 

a = e (25) 

R{a) = bo + b{a) . 

Further, by considering equation (28), we see that the curvature of this shape will always 
be positive for sufficiently large bo independent of the form of b{6). 
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V. CONCLUSIONS 


We have demonstrated that every two-dimensional rigid body making elastic collisions 
between parallel walls is equivalent to a traditional billiard problem in which a point parti¬ 
cle makes specular reflections with walls of uniquely determined shape. Certain shapes of 
particular geometrical interest may encourage the consideration of specific billiard problems 
which might otherwise have escaped attention. Rectangles and shapes of constant width 
may warrant investigation. Strongly asymmetric rigid bodies strike us as being of particular 
interest. We note that the present results can be extended to include the effects of external 
fields. 

Similarly, every such billiard problem can be associated with a family of equivalent rigid 
body problems. This offers the possibility of viewing the extensive billiard literature in a new 
and potentially revealing light. This equivalence can be extended. It is evident, for example, 
that similar results can be obtained for a two-dimensional body confined in a rectangular 
box. There is, of course, nothing to preclude the consideration of billiard problems in higher 
dimensional spaces. We believe that the existence of simple physical systems can give such 
systems an intuitive immediacy which they might otherwise lack. Also, the motion of a rigid 
body in three and higher dimensions is associated with a non-commutative group. Thus, we 
expect it to be equivalent to a new type of billiard motion. 

Finally, the analogous quantum mechanical problem of rigid bodies confined between 
parallel walls is of evident interest and is a current focus of our attention. Here, the features 
of the quantal results require a new interpretation, due to the new association with a different 
physical property, the tunneling from heads to tails. 

We would like to acknowledge helpful discussions with Dr. Serge Troubetzkoy. One of 
us (R.C.) would like to acknowledge helpful discussions with M. Khanji. This work was 
supported in part by the US Department of Energy under grant No. DE-FG02-88ER 40388 
(A.D.J. and R.C.) and by the National Science Foundation (N.L.B.). 
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VI. APPENDIX 


Here we shall start from the wall shape obtained for an elliptic disc and consider the 
various rigid body shapes which result from different choices of h^. Specihcally, we take 

■l + ecos2 0^^/^ 


h{e) = &o + 


1+6 


( 26 ) 


In the vicinity of a = 0, we hnd 


R{a) 


ipQ + 1 ) 


1 (1 + 6o)e 2 

2 l + 6o(l + e) " 


(27) 
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It is useful to recall that the curvature of R{a) is 


+ 2R'^ - RR" 

{R^ + /?' 2 ) 3/2 


(28) 


where primes indicate derivatives with respect to a. The curvature corresponding to equation 
(24) at a = 0 is thus 


A'(0) 


1 + 6 

1 + 6o(l + e) 


(29) 


Evidently, there is a critical value of 6o, be = —1/(1 + e) such that -ft"(0) is negative for 
bo < be- Since the wall function was assumed to come from the convex hull of the original 
rigid body, this result is unphysical. There is no corresponding rigid body for bo < be- It is 
instructive to describe what happens as one integrates from Q; = 7r/2toQ; = 0in order to 
obtain R{a): 

Case I. ( b > be) Here, R{a) approaches bo + 1 quadratically (as expected) and R'{a) 
approaches 0 linearly as a —> 0. The result is a physical shape which depends on bo. The 
result is consistent in the sense that the original wall function can be reconstructed from 
this R{a) using the general prescription of section II. 

Case II. ( b = be) This is the limiting physical case. Now, R{a) approaches bo + I like 
and R'{a) goes like Clearly, R" is badly behaved at a = 0, but this causes no 

difficulty. Again, the result is a completely physical shape which permits reconstruction the 
original wall function. This critical shape is shown in Fig.4. 

Case III. ( b < be) R{a) approaches 6o + 1 linearly, but R'{a) does not approach 0. 
Rather, it approaches some non-zero value, R'q, linearly. The putative shape has a discon¬ 
tinuous derivative at a = 0 (and a = tt). To see what this means, return to equation (12). 
The range 0+ < a < 7r/2 corresponds to covering the range Oc < 0 < n/2 with 


9c = tan ^ 


R'o \ 
(J'o + l)/ 


(30) 


The hnite region between —9e < 6 < +6e is all governed by the inhnitesimal region 0_ < 
CK < 0+. There is nothing wrong with a single point on the body determining the wall shape 
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for a finite angular interval. This occurs whenever the rigid body has a corner and R{a) has 
a discontinuous derivative. However, as in equation (5), this should correspond to a shape 
{bo + 1)1 cos0|. Evidently, this is not the shape of bo + b{6) chosen here nor will it be for 
other choices of b{9). Hence, the resulting body shape will fail to reproduce bo + b{6) over 
the range —9(. < 9 < 9c, and thus R{a) is not acceptable. We remark that, in the present 
example, bo + b{9) remains positive (and thus not manifestly unphysical) for values of bo as 
small as —1/Vl + e which is less than the limiting value of be quoted above. 

Similar estimates on be based on the existence of regions of negative curvature at the 
extrema of b{9) can be made for arbitrary (differentiable) wall functions since these extrema 
necessarily coincide with the extrema of R{a). Negative curvature always indicates an 
unphysical shape. While we have not demonstrated that the hrst appearance of negative 
curvature will be at an extremum, these ease of these estimates suggests that they may be 
useful. 
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FIGURES 

FIG. 1. Phase portrait (in units of tt) of the contact point 6c versus the angle of incidence (/>, 


with respect to the local normal 9c for the elliptic system e = 0.03 and H = 3. 

FIG. 2. Phase space diagram for the billiard system with e = 0.55 and H = 1.9. 

FIG. 3. The rigid body shape {R{a)) corresponding to the limiting physical 
bo = bc = -1/2, (e = 1). 


case of 
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